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INTRODUCTION 


Electromechanical  transduction  taking  place  via  the  piezoelectric 
effect  is  characterized  phenomenologically  by  constitutive  equa¬ 
tions  that  relate  the  elastic  and  electric  variables.  These  equa¬ 
tions  take  a  variety  of  forms,  depending  upon  the  choice  of  inde¬ 
pendent  and  dependent  variables;  the  choice  is  normally  dictated  by 
the  application.  For  example,  piezoelectric  resonators  in  the  form 
of  thickness  mode  plates  are  most  easily  treated  using  the  isagric 
elastic  stiffnesses  [cE],  the  piezoelectric  stress  constants  [e] , 
and  the  dielectric  permittivities  at  constant  strain  [(eps)S]. 

Various  measurement  techniques  yield  values  for  the  elements  of  a 
particular  coefficient  set  more  directly  than  those  of  another. 
The  coefficients  appearing  in  the  different  equation  sets  are, 
however,  interrelated,  so  that  once  any  one  complete  set  is  avail¬ 
able,  all  the  other  sets  of  elements  may  be  found.  The  most  accu¬ 
rate  and  precise  experimental  results  to  date  have  been  from  plate 
resonator  (resonance)  and  pulse-echo  (transit-time)  measurements. 
From  the  [cE],  [e],  and  [(eps)S]  matrices  determined  therefrom, 
those  matrices  representing  material  properties  expressed  in  the 
other  alternative  forms  may  be  calculated. 

Electrooptical  applications  are  becoming  increasingly  important. 
So  also  are  treatments  of  piezoelectric  and  ferroelectric  phenomena 
from  the  standpoint  of  molecular  interactions.  In  both  of  these 
cases  the  constitutive  equations  using  polarization  as  the  inde¬ 
pendent  electrical  variable,  rather  than  either  electric  intensity 
or  displacement,  assume  greater  importance  than  the  sets  tradition¬ 
ally  used  for  transducer,  signal  processing,  and  resonator  applica¬ 
tions. 

In  this  report  we  give  the  complete  sets  of  linear  constitutive 
equations  relating  elastic  and  electric  fields.  For  each  equation 
set  the  numerical  values  are  computed  for  lithium  niobate,  the 
preeminent  electrooptic  and  acoustooptic  material,  from  the  meas¬ 
ured  [cE] ,  [e],  and  [(eps)S]  values  of  Smith  and  Welsh  (Ref.l). 
Coupling  to  the  thermal  field  is  neglected.  Rationalized  mks  units 
are  used  throughout. 


CONSTITUTIVE  EQUATION  SETS 

Symbols  and  units  for  the  quantities  employed  are  given  in  Table  1. 
In  terms  of  these,  six  constitutive  equation  sets  are  used.  Of 
these,  electric  intensity,  dielectric  displacement,  and  polariza¬ 
tion  each  appear  in  two  sets  as  an  independent  variable.  The  sets 
are,  in  compressed  matrix  notation,  as  follows.  A  prime  denotes 
transpose;  [I]  is  the  unit  matrix. 


I.  The  Piezoelectric  Stress  Constant  Set 


(  1) 
(  2) 
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[T]  =  (CE]  [S]  -  (e]'  [E] 

(D)  =  [  e)  [S]  +  [ (eps)S]  [E] 


TABLE  1.  SYMBOLS,  UNITS,  AND 

DEFINITIONS. 

QUANTITY 

UNIT 

SYMBOL/DEFINITION 

Elastic  stress 

N/m2 

[T] 

Elastic  strain 

— 

[S] 

Electric  intensity 

V/m 

[E] 

Dielectric  displacement 

C/m2 

[D] 

Dielectric  polarization 

C/m2 

[P] 

Elastic  compliance  at 
constant  [E],  [D],  [P] 

m2/N 

[CE],  [CD],  [CP] 

Elastic  stiffness  at 
constant  [E] ,  [D],  [P] 

N/m2 

[sE],  [sD],  [sP] 

Dielectric  permittivity 
at  constant  [T],  [S] 

F/m 

[ (eps) T] ,  [ (eps) S] 

Dielectric  constant, 
relative,  at 
constant  [T],  [S] 

[(  Kr ) T ] ,  [(  Kr) S] 
=[ (eps)T]/ (eps) o, 

[ (eps) S]/ (eps) o 

Dielectric  impermeability 
at  constant  [T] ,  [S] 

m/F 

[(bet)T],  [fbet)S] 
-[ (eps)T]  -f  , 
t (eps) S]  ( 

Dielectric  impermeability, 
relative,  at  constant 
[T],  [S] 

[ (betr) T] ,  [ (betr)S] 
=[ (bet)T]*(eps)o, 

[ (bet) S] * (eps) o 
=[(  Kr) T]  , 

[ (  Kr) S]  l-1* 

Dielectric  susceptibility 
at  constant  [T] ,  [S] 

F/m 

[ (chi) T] ,  [ (chi)S] 
=[(  Kr)T-I) * (eps) o, 

[(  Kr)S-I]*(eps)o 

Dielectric  susceptibility, 
relative,  at  constant 
[T],  [S] 

[ (chir)T) ,  [ (chir) S] 
=[ (chi)T]/(eps)o, 

[ (chi) S]/ (eps) o 

Reciprocal  dielectric 
susceptibility  at 
constant  [T] ,  [S] 

m/F 

[ (zet)TJ,  [(zet)S] 

= [ (chi) T] 

[ (chi) S]  <-1' 

Reciprocal  dielectric 
susceptibility,  relative, 
at  constant  [T] ,  [S] 

[(zetr)T],  [(zetr)S) 
=[ (zet)T]*(eps)o, 

[  (zet)S]*(eps)o 

Piezoelectric  stress 
constant 

C/m2 

2 

te] 

TABLE  1.  SYMBOLS,  UNITS,  AND  DEFINITIONS,  (continued) 


QUANTITY 


UNIT  SYMBOL/DEFINITION 


Piezoelectric 

coefficient 

strain 

m/V  =  C/N 

[d] 

Piezoelectric 

modulus 

stress 

N/C  =  V/m 

[h] 

Piezoelectric 

constant 

strain 

m2/C 

[g] 

Piezoelectric 

modulus 

polarization 

V/m  =  N/C 

[a] 

Piezoelectric 

constant 

polarization 

m2/C 

[b] 

Note:  Square  brackets,  sic:  [  ],  denote  matrices. 


3 


II.  The  Piezoelectric  Strain  Coefficient  Set 

[S]  =  [SE]  [T]  +  [d]'  [E]  (  3) 

[D]  =  [  d]  [T]  +  [ (eps)T]  [E]  (  4) 

III.  The  Piezoelectric  Stress  Modulus  Set 

[T]  =  [CD]  [S]  -  [h]  >  [D]  (  5) 

[E]  =  -[  h]  [S]  +  [  (bet) S]  [DJ  (  6) 

IV.  The  Piezoelectric  Strain  Constant  Set 

[S]  =  [sD]  [T]  +  Cg] #  [D]  (  7) 

[E]  =  -[  g]  [T]  +  [ (bet)T]  [D]  (  8) 

V.  The  Piezoelectric  Polarization  Modulus  Set 

[T]  =  [CP]  [S]  -  [a]'  [P]  (  9) 

[E]  =  -[  a]  [S]  +  [  (zet)S]  [P]  (10) 

VI.  The  Piezoelectric  Polarization  Constant  Set 

[S]  *  [sP]  [T]  +  [b]'  [P]  (11) 

[E]  =  -[  b]  [T]  +  [  (zet)T]  [P]  (12) 

The  electric  variables  are  connected  by  the  relation 

[D]  =  (eps)o  *  [E]  +  [P]  (13) 

where  (eps)o  is  the  permittivity  of  free  space,  defined  by 

(eps)o  *  (mu)o  *  (c)  *  (c)  =  1  ;  (14) 


(mu)o  is  .the  permeability  of  free  space,  equal,  by  definition,  to 
4  PI  *  10(-7>,  and  (c)  is  the  velocity  of  light  in  vacuo  and,  also 
by  definition,  is  equal  exactly  to  2.99792458  x  108  m/s. 

From  (13)  the  expressions  for  the  remaining  electric  variables 
associated,  respectively,  with  the  six  equation  sets  (  1)  to  (12) 
may  be  found: 


[P]  = 

[  e]  [S] 

+  [ (chi) S] 

[E] 

(15) 

CP]  - 

[  d]  [T] 

+  [ (chi)T] 

[E] 

(16) 

[P]  - 

(eps)o  * 

[  h)  [S)  + 

[I  - 

(eps)o  * 

(bet)S] 

[D] 

(17) 

[P]  = 

(eps)o  * 

t  g]  IT]  + 

[I  - 

(eps)o  * 

(bet)T] 

[D] 

(18) 
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[D] 

=  - (eps) o  * 

[  a] 

[S]  + 

[I  + 

(eps)o  * 

(zet)S] 

[P] 

(19) 

CD] 

=  - (eps) o  * 

C  b] 

[T]  + 

[I  + 

(eps)o  * 

(zet) T] 

[P] 

(20) 

RELATIONS  AMONG  MATERIAL  CONSTANTS 

The  material  constants  are  interrelated  by  the  following  formulas: 

[cX]  [sX]  =  [ (eps) Y]  [ (bet) Y]  =  [I]  (21) 

[ (chi) Y)  [ (zet) Y]  =  [(  Kr) Y  -  (chir)Y]  =  [I]  (22) 

In  (21)  and  (22),  X  =  E,  D,  or  P  and  Y  =  T  or  S. 

(cD]  -  [cE]  =  Ch] •  [e]  =  [e]'  [hj 

-  [h]'  [ (eps) S]  [h]  =  [e] '  [ (bet) S]  [e] 

=  [a]'  [e  -  h  *  (eps)o]  =  [e  -  h  *  (eps)o]'  [a]  (23) 

[cP]  -  [cD]  =  [h] '  [a]  *  (eps)o  *  [a]'  [h]  *  (eps)o 

=  [h] *  [ (eps) S]  [ (zet) S]  [h]  *  (eps) o 

=  [a]*  [ (bet) S]  [ (chi) S]  [a]  *  (eps)o 

=  [a  -  h] '  fe]  =  re]'  [a  -  h]  (24) 

[cP]  -  [cE]  =  [a]'  [e]  «  [e]'  [a] 

=  [a]'  [ (chi) S]  [a]  =  [e] •  [(zet)S]  [e] 

=  [h] '  [e  +  a  *  (eps)o]  =  [e  +  a  *  (eps)o]'  [h]  (25) 

[sE]  -  [sD]  =  [d]»  [g]  =  [g]'  [d] 

-  [d] •  [ (bet) t]  [d]  =  [g] ■  [ (eps) t]  [g] 

=  tb] '  [d  -  g  *  (eps)o]  =  [d  -  g  *  (eps)o]'  [b]  (26) 

[sD]  -  [sP]  =  [b] •  [g]  *  (eps) o  =  [g] •  [b]  *  (eps)o 

=  [g]‘  [ (eps) T]  [ (zet)T]  [g]  *  (eps)o 

=  [b] ■  [ (bet) T]  [ (chi) T]  [b]  *  (eps) o 

=  [b  -  g] •  [d]  =  [d]»  [b  -  g]  (27) 

[sE]  -  [sP]  =  [b] •  [d]  *  [d] •  [b] 

=  [b] •  [ (chi) T]  [b]  =  [d]«  [ (zet)T]  [d] 

=  [g] •  [d  +  b  *  (eps)o]  =  [d  +  b  *  (eps)o]'  [g]  (28) 
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t (zet)S]  -  ( (zet) T]  =  [b]  [a]'  =  [a]  [b] • 

=  [b]  [cP]  [b] •  =  [a]  [sP]  [a]*  (29) 

[ (chi) T]  -  [ (chi) S]  -  [(eps)T]  -  [<eps)S] 

=  [e]  [d]'  =  [d]  [e]' 

«  [d]  [cE)  [d] •  =  [e]  [sE]  [e] '  (30) 

[ (bet) S]  -  [ (bet)T]  =  [h]  [g] •  =  [g]  [h] • 

-  [g]  (CD]  [g]'  =  [h]  (sD)  [h] '  (31) 

[e]  =  [d]  [cE]  =  [ (eps) S]  [h]  =  [(chi)S]  [a]  (32) 

[d]  =  [e]  [sE]  =  [ (eps) T]  [g]  =  [(chi)T]  [b]  (33) 

[h]  =  [g]  [CD]  =  [ (bet) S]  [e]  =  [(chi)S]  [(bet)S)  [a] 

=  [I  -  (bet)S  *  (eps)o]  [a]  (34) 

[g]  =  [h]  [sD]  =  [ (bet) T]  [d]  =  [(chi)T]  [(bet)T]  [b] 

=  [I  -  (bet) T  *  (eps) o]  [b]  (35) 

[a]  =  [b]  [cP]  =  [ (zet) S]  [e]  =  ((eps)S]  [(zet)S]  [h] 

=  [I  +  ( zet) S  *  (eps)o]  [h]  (36) 

[b]  =  [a]  [sP]  =  [ { zet) T]  [d]  =  [(eps)T)  [(zet)T]  [g] 

=  [I  +  (zet)T  *  (eps) o ]  [g]  (37) 


Some  alternative  relations  are  the  following: 


[a  -  h]  =  [ (zet)S]  [h]  *  (eps)o 

«  [ (bet) S]  [a]  *  (eps)o  (38) 

(b  -  g]  =  [ (zet) T]  [g]  *  (eps)o 

=  [ (bet)T]  [b]  *  (eps) o  (39) 

[e  +  a  *  (eps) o]  =  [(eps)S]  [a]  (40) 

[d  +  b  *  (eps)o]  =  ((eps)T]  [b]  (41) 

[e  -  h  *  (eps)o]  =  [(chi)S]  [h]  (42) 

[d  -  g  *  (eps)o]  =  [(chi)T]  [gj  (43) 


Equations  (21)  to  (43)  result  from  equating  like  dependent  vari 
ables  pairs  selected  from  equations  (  1)  to  (12)  and  (15)  to  (20) 
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Each  pair  yields  one  equation  in  three  variables,  one  mechanical 
and  two  electrical,  or  vice  versa.  Two  other  equations  exist, 
again  from  (1)  to  (12)  and  (15)  to  (20),  that  contain  the  same 
three  variables  found  in  each  paired  equation.  One  of  these  auxil¬ 
iary  equations  is  used  to  eliminate  one  of  the  two  variables  of  the 
same  kind;  the  result  is  one  equation  in  two  variables,  one  elec¬ 
trical  and  one  mechanical.  These  are  now  independent  variables,  so 
the  coefficients  must  vanish;  two  relations  between  the  material 
coefficients  result.  As  an  example,  (  3)  and  (  7)  both  have  [S]  as 
dependent  variable.  Equating  them  produces  one  relation  in  [T] , 
[E],  and  (D] ;  one  of  the  electrical  variables  must  be  eliminated. 
This  is  done  by  using  either  (  4)  or  (  8) ;  each  contains  the  same 
three  variables.  If  (  8)  is  used  to  eliminate  [E],  one  obtains  [sE 
-  d'  g  -  sD]  [T]  =  [d*  (bet)T  -  g1]  [D] .  Therefore,  [sE]  -  [sD]  = 
[d] '  [g]  and  [g]  =  [(bet)T]  [d] .  Use  of  (  4)  instead  of  (  8)  leads 
to  the  equations  [sE]  -  [sD]  =  [g] '  [d]  and  [d]  =  [(eps)T]  [g]. 
There  are  36  pairs,  six  each  equating  [S]  and  [T] ,  and  eight  each 
equating  [ E ] ,  [D],  and  [ P ] .  The  72  relations  contain  many  redun¬ 
dancies.  Relations  between  the  elastic,  piezoelectric,  and  dielec¬ 
tric  constants  are  shown  schematically  in  Tables  2  and  3. 


CALCULATION  SEQUENCE 

Using  as  input  [cE],  [e],  and  [(eps)S],  one  may  compute  the  remain¬ 
ing  quantities  in  a  variety  of  ways.  The  following  sequence  is 


typical: 

[sE]  =  [cE]  (~1)  (44) 
t (bet) S]  =  [ (eps) S]  f’1)  (45) 
[d]  =  [e]  [sE]  (46) 
[h]  =  ( (bet) S]  [e]  (47) 
t (eps)T]  -  [ (eps) S]  =  [e]  [d] •  (48) 
C (eps) T]  =  [ (eps) S]  +  [e]  [d]  '  (49) 
[ (bet) T]  =  [ (eps) T]  (50) 
[CD]  -  [cE]  =  [e]'  [h]  (51) 
[CD]  =  [cE]  +  [e]'  [h]  (52) 
[g]  =  [ (bet) T]  [d]  (53) 
[sE]  -  [SD]  =  [d] '  [g]  (54) 
[sD]  =  [sE]  -  [d] '  [g]  (55) 
[ (betr) S]  =  [ (bet) S]  *  (eps)o  (56) 
[ (zetr)S]  =  [ (betr)S]  (I  -  (betr)S]  f"1)  (57) 
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t (zet)S]  =  [(zetr)]  /  (eps)o  (58) 
[(betr)T]  =  [(bet)T]  *  (eps)o  (59) 
( (zetr)T]  *  [ (betrT]  [I  -  (betr)T]  (60) 
[(zet)T  =  [(zetr)T]  /  (eps)o  (61) 
( (chi) S]  *  [  (zet) S]  (-1>  (62) 
[  (chi)T]  =  [ (zet)T]  (63) 

(a]  *  [ (zet)S]  [e]  (64) 

[b]  *  [ (zet)T]  [d]  (65) 
(cP]  -  [CE]  =  [e]'  [a)  (66) 
[cP]  =  [cE]  +  [e]'  [a]  (67) 
[cP]  -  [cD]  »  [a]'  [h]  *  (eps)o  (68) 
[sE]  -  [sP]  =  [d]  •  [b]  (69) 
[sP]  =  [sE]  -  [d] *  [b]  (70) 
[SD]  -  [sP]  =  [g]'  [b]  *  (eps)o  (71) 
[ (bet) S]  -  [ (bet)T]  =  [h]  [g] •  (72) 
[ (chi)T]  -  [ (chi) S]  -  [ (eps)T]  -  [(eps)S]  (73) 
[  (zet) S]  -  [  (zet) T]  =  [a]  [b]  '  (74) 


A  number  of  these  relations  are  used  as  checks.  For  example, 
[(bet)S]  and  [(bet)T]  are  known  from  (45)  and  (50),  but  the  differ¬ 
ence  is  recomputed  in  (72) . 

EXPLICIT  FORMULAS  FOR  POINT  GROUP  3m 
Elastic: 

The  6x6  elastic  constant  portions  of  Tables  4  and  5  partition  into 
4x4  and  2x2  submatrices.  The  4x4  elastic  stiffness  and  compliance 
submatrices  are  interrelated  by  formulas  (75)  to  (93) ,  taken  from 
Cady  (Ref.  2): 


A 

2 

S33 

* 

(Sll  +  Sl2 )  - 

2 

* 

Sl3 

*  Sl3 

(75) 

B 

2 

s44 

* 

(Sll  -  sl2 )  - 

2 

* 

s!4 

*  Sl4 

(76) 

2 

* 

Cll 

2 

s33  /  A  +  s44 

/ 

B 

(77) 

2 

* 

Cl2 

2 

s33  /A  -  s44 

/ 

B 

(78) 
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TABLE  4.  ELASTOPIEZODIELECTRIC  MATRICES  FOR  POINT  GROUP  3m: 


THE 

m 

0) 

t_j 

[h]. 

AND  [a] 

SETS. 

11 

12 

13 

14 

00 

00 

] 

00 

-22 

31 

3 

cE 

3  e ' 

12 

11 

13 

-14 

00 

00 

3 

00 

22 

31 

—  -- 

■3 - 

3 

e 

3 (eps)S 

13 

13 

33 

00 

00 

00 

3 

1 

00 

00 

33 

14 

-14 

00 

44 

00 

00 

J 

3 

00 

15 

00 

3 

cD 

3  h' 

00 

00 

00 

00 

44 

14 

3 

15 

00 

00 

— 

3 - 

3 

h 

] (bet) S 

00 

00 

00 

00 

14 

66 

3 

-22 

00 

00 

00 

00 

00 

00 

15 

-22 

J 

3 

11 

00 

00 

3 

CP 

3  a* 

-22 

22 

00 

15 

00 

00 

3 

00 

11 

00 

— 

3 - 

3 

a 

] (zet) S 

31 

31 

33 

00 

00 

00 

3 

00 

00 

33 

66 

= 

u 

l  II 
II 

rH  || 
H  || 

12)  / 

2 

Matrix  entries  show  only  subscripts. 


TABLE  5.  ELASTOPIEZODIELECTRIC  MATRICES  FOR  POINT  GROUP  3m: 


THE 

[d], 

[g]. 

AND  [b] 

SETS. 

11 

12 

13 

14 

00 

00  ] 

00 

-22 

31 

3 

SE 

3  d' 

12 

11 

13 

-14 

00 

00  ] 

00 

22 

31 

— 

3 - 

3 

d 

] (eps)T 

13 

13 

33 

00 

00 

00  ] 

00 

00 

33 

14 

-14 

00 

44 

00 

J 

00  ] 

00 

15 

00 

3 

SD 

3  g' 

00 

00 

00 

00 

44 

14*2  ) 

15 

00 

00 

— 

•3 - 

3 

g 

] (bet) T 

00 

00 

00 

00 

14*2  66  ] 

1 

-22*2 

00 

00 

00 

00 

00 

00 

15 

j 

-22*2  ] 

11 

00 

00 

3 

SP 

3  *>• 

-22 

22 

00 

15 

00 

00  ] 

00 

11 

00 

— 

■3 - 

3 

b 

] (zet)T 

31 

31 

33 

00 

00 

00  ) 

00 

00 

33 

66  =  (11  -  12)  *  2 

Matrix  elements  show  only  subscripts. 
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(79a),  (79b) 


cl’  =  -  Sl3  /  A  ;  cl4  =  -  sl4  /  B 

C33  =  (sll  +  Sl2)  /  A 

C44  =  (sll  -  Sl2 )  /  B 

C66  =  (Cll  -  Cl2)  /  2  =  S44  /  (2  *  B) 

K  =  C33  *  (cll  +  Cl2 )  -  2  *  cl3  *  Cl3 

L  =  C44  *  (cll  -  Cl2 )  -  2  *  Cl4  *  Cl4 

2  *  sll  =  c33  /  K  +  C44  /  L 

2  *  sl2  ■  c33  /  K  -  c44  /  L 

sl3  =  -  Cl3  /  K  ;  Sl4  =  -  Cl4  /  L 

S33  =  (cll  +  Cl2 )  /  K 

S44  =  (cll  -  Cl2 )  /  L 

S66  =  (sll  -  sl2 )  =  2  *  C44  /  L 

det  (4x4)  [s]  «  A  *  B 

det  (4X4)  [c]  =  K  *  L 

A*K  =  B*L  =  A*B*K*L=1 


(87a) ,  (87b) 


(79) 

(80) 
(81) 
(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 

(89) 

(90) 

(91) 

(92) 

(93) 


Formulas  (75)  to  (93)  hold  for  each  set  of  constant  electrical 


conditions:  either  E,  D,  or  P  constant. 

[cD]  -  [cE]  =  [del  cDE]  =  [e]'  [h]  =  [h] '  [e]  (23) 

del  cDEll  =  +  e22  h22  +  e31  h31  (94) 

del  CDE12  =  -  e22  h22  +  e31  h31  (95) 

del  CDE13  =  +  e31  h33  =  +  h31  e33  (96) 

del  CDE14  =  -  e22  hl5  =  -  h22  el5  (97) 

del  CDE33  =  +  e33  h33  (98) 

del  CDE44  =  +  el5  hl5  (99) 

del  CDE66  =  +  e22  h22  (100) 

[cP]  -  [cD]  =  [del  cPD]  =  [a] •  [h]  *  (eps)o 

=  [h] *  [a]  *  (eps)o  (24) 
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del  cPDll  -  (  +  a22  h22  +  a31  h31  )  *  (eps)o  (101) 

del  cPD12  =  (  -  a22  h22  +  a31  h31  )  *  (eps)o  (102) 

del  cPD13  =  (  +  a31  h33  )  *  (eps)o 

=  (  +  h31  a33  )  *  (eps)o  (103) 

del  cPD14  =  (  -  a22  hl5  )  *  (eps)o 

=  (  -  h22  al5  )  *  (eps)o  (104) 

del  cPD33  =  (  +  a33  h33  )  *  (eps)o  (105) 

del  cPD44  =  (  +  al5  hl5  )  *  (eps)o  (106) 

del  CPD66  =  (  +  a22  h22  )  *  (eps)o  (107) 

[cP]  -  [cE]  -  [del  cPE]  =  [e]  •  [a)  =  [a]'  [e]  (25) 

del  cPEll  =  +  e22  a22  +  e31  a31  (108) 

del  cPE12  =  -  e22  a22  +  e31  a31  (109) 

del  CPE13  -  +  e31  a33  =  +  a31  e33  (110) 

del  CPE14  =  -  e22  al5  =  -  a22  el5  (ill) 

del  CPE33  =  +  e33  a33  (112) 

del  cPE44  =  +  el5  al5  (113) 

del  CPE66  =  +  e22  a22  (114) 

From  the  del  cl3  entries  we  have  the  ratios 

e31  /  e33  =  h31  /  h33  =  a31  /  a33.  (115) 

From  the  del  cl4  entries  we  have  the  further  ratios 

el5  /  e22  =  hl5  /  h22  =  al5  /  a22.  (116) 

[sE]  -  (sD]  =  [del  sED]  =  [d]  •  [g]  =  [g] '  [d]  (26) 

del  sEDll  =  +  d22  g22  +  d31  g31  (117) 

del  SED12  =  -  d22  g22  +  d31  g31  (118) 

del  SED13  =  +  d31  g33  =  +  g31  d33  (119) 

del  sED14  =  -  d22  gl5  =  -  g22  dl5  (120) 

del  SED33  =  +  d33  g33  (121) 
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del  SED44  =  +  dl5  gl5  (122) 

del  SED66  =  +  d22  g22  *  4  (123) 

[sD]  -  [sP)  =  [g]'  [b]  *  (eps)o 

=  [b] •  [g]  *  (eps)o  (27) 

del  sDPll  =  (  +  g22  b22  +  g31  b31  )  *  (eps)o  (124) 

del  SDP12  =  (  -  g22  b22  +  g31  b31  )  *  (eps)o  (125) 

del  sDP13  =  (  +  g31  b33  )  *  (eps)o 

=  (  +  b31  g33  )  *  (eps)o  (126) 

del  SDP14  =  (  -  g22  bl5  )  *  (eps)o 

=  (  -  b22  gl5  )  *  (eps)o  (127) 

del  sDP33  =  (  +  g33  b33  )  *  (eps)o  (128) 

del  sDP44  =  (  +  gl5  bl5  )  *  (eps)o  (129) 

del  SDP66  =  (  +  g22  b22  )  *  4  *  (eps)o  (130) 

(sE]  -  [sP]  =  [del  sEP)  =  [b] '  [d]  =  [d] '  [b]  (28) 

del  sEPll  *  +  d22  b22  +  d31  b31  (131) 

del  SEP12  =  -  d22  b22  +  d31  b31  (132) 

del  SEP13  -  +  d31  b33  -  +  b31  d33  (133) 

del  SEP14  =  -  d22  bl5  =  -  b22  dl5  (134) 

del  SEP33  =  +  d33  b33  (135) 

del  sEP44  =  +  dl5  bl5  (136) 

del  SEP66  =  +  d22  b22  *  4  (137) 

From  the  del  sl3  entries  we  have  the  ratios 

d31  /  d33  =  g31  /  g33  =  b31  /  b33.  (138) 

From  the  del  sl4  entries  we  have  the  further  ratios 

dl5  /  d22  =  gl5  /  g22  =  bl5  /  b22.  (139) 

Piezoelectric: 

[d]  *  [e)  [sE]  (33) 
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dl5  *  +  el5  sE44  -  e22  sE14  *  2  (140) 
d22  -  +  e22  (sEll  -  sE12)  -  el5  sE14  (141) 
d31  -  +  e31  (sEll  +  sE12 )  +  e33  sE13  (142) 
d33  =  +  e33  sE33  +  el3  SE13  *  2  (143) 
[h]  =  [ (bet) S)  [e]  (34) 
hl5  =  (bet)Sll  el5  (144) 
h22  =  (bet)Sll  e22  (145) 
h31  =  (bet) S33  e31  (146) 
h33  -  (bet) S33  e33  (147) 
[g]  =  [ (bet) T]  [d]  (35) 
gl5  =  (bet) Til  dl5  (148) 
g22  =  (bet) Til  d22  (149) 
g31  =  (bet) T33  d31  (150) 
g33  =  (bet) T33  d33  (151) 

[a]  =  [ (zet)S]  (e]  (36) 
al5  =  (zet)Sll  el5  (152) 
a22  =  (zet)Sll  e22  (153) 
a31  =  (zet)S33  e31  (154) 
a33  =  (zet)S33  e31  (155) 

[b]  =  [ (zet)T]  (d]  (37) 
bl5  =  (zet)Tll  dl5  (156) 
b22  =  (zet)Tll  d22  (157) 
b31  =  (zet)T33  d31  (158) 
b33  -  (zet)T33  d33  (159) 

Dielectric: 

[ (bet) Y)  =  [(eps)Y](_1)  (21) 
(bet)Yll  *  1  /  (eps)Yll  (16°) 
(bet)  Y33  *  1  /  (eps)  Y33  (161) 
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[ (zetr) Y]  =  [ (betr) Y]  [I  -  (betr)Y]<_1>  (162) 

(zet)Yll  =  1  /  ((eps)Yll  -  (eps)o)  (163) 

(zet)Y33  =  1  /  ( (eps) Y33  -  (eps)o)  (164) 

[(eps)T  -  (eps)S]  =  [del  (eps)]  =  [e]  [d] •  = 

[ (chi)T  -  (chi) S]  =  [del  (chi)]  =  [d]  [e] •  (30) 

del  (eps) 11  =  del  (chi) 11  =  +  el5  dl5  +  e22  d22  *  2  (165) 

del  (eps) 33  =  del  (chi) 33  =  +  e33  d33  +  e31  d31  *  2  (166) 

[  (bet) S  -  (bet) T]  =  [h]  [g]  •  =  [g]  [h]  •  (31) 

del  (bet) 11  =  +  hl5  gl5  +  h22  g22  *  2  (167) 

del  (bet) 33  =  +  h33  g33  +  h31  g3l  *  2  (168) 

[ (zet) S  -  ( zet) T]  =  [del  (zet)]  =  [a]  [b] •  =  [b]  [a]'  (169) 

del  (zet) 11  =  +  al5  bl5  +  a22  b22  *  2  (170) 

del  (zet) 33  =  +  a33  b33  +  a31  b31  *  2  (171) 
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PIEZOELECTRIC  CONSTANT  DETERMINATIONS 


Determinations  of  the  [e] ,  [d] ,  [h],  and  [g]  values  for  lithium 
niobate  have  been  made  by  various  investigators.  The  most  consist¬ 
ent  of  these  are  shown  in  Tables  6,  7,  8,  and  9  for  comparison 
purposes . 


TABLE  6.  PIEZOELECTRIC  STRESS  CONSTANT  [e] . 


el5 

e33 

T  (OC) 

Ref 

3.7 

2.5 

0.2 

1.3 

RT 

(3) 

3.60 

2.52 

0.747 

1.67 

RT 

(4) 

3.76 

2.43 

0.23 

1.33 

25 

(1) 

3.8 

2.5 

0.35 

1.42 

20 

(5) 

3.83 

2.37 

1.80 

RT 

(6) 

3.61 

2.40 

0.28 

1.59 

RT 

(7) 

Estimated 

inaccuracies:  el5, 

e22 :  3% 

;  e31:  60% 

;  e33 :  15%. 

Units:  C/m2. 


TABLE 

7.  PIEZOELECTRIC  STRAIN 

COEFFICIENT 

[d]. 

dl5 

d22 

d31 

d33 

T  (OC) 

Ref 

74.0 

20.8 

-0.863 

16.2 

20 

(8) 

68. 

21. 

-1. 

6. 

RT 

(3) 

64.3 

20.6 

+  1.15 

6.53 

RT 

(4) 

69.2 

20.8 

-0.85 

6.0 

25 

(1) 

68.9 

20.9 

+0.003 

5.8 

20 

(5) 

65.36 

20.29 

-1.22 

8.27 

RT 

(7) 

Units: 

10(-12>  C/N. 

TABLE  8.  PIEZOELECTRIC  STRESS  MODULUS  [h] . 


hl5 

h22 

h33 

T  (OC) 

Ref 

9.5 

6.4 

0.8 

5.1 

RT 

(3) 

9.04 

6.33 

3.07 

6.86 

RT 

(4) 

9.59 

6.20 

0.93 

5.38 

25 

(1) 

9.2 

6.1 

1.3 

5.87 

20 

(5) 

9.16 

6.09 

1.2 

6.88 

RT 

(7) 

Units:  10 (  9>  N/C 
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T  (OC) 


Ref 


TABLE  9. 
gl5 


PIEZOELECTRIC  STRAIN  CONSTANT  [g] . 


g22 


g3l 


g33 


91. 

28. 

-4. 

87.7 

28.1 

+4.79 

91.8 

27.6 

-3.3 

92.5 

28.1 

+0.01 

89.5 

27.8 

-4.87 

Units: 

10("  3)  m2/C. 

23. 

RT 

(3) 

24.4 

RT 

(4) 

23.6 

25 

(1) 

23.5 

20 

(5) 

33.0 

RT 

(7) 

INPUT  VALUES  FOR  LI  NB  03 

The  values  measured  by  Smith  and  Welsh  (Ref.  1)  using  the  pulse 
echo  (transit-time)  technique  are  as  follows: 


TABLE  10.  ISAGRIC  ELASTIC  STIFFNESSES. 

CE11  CE12  CE13  CE14  CE33  CE44  CE66 

203.0  57.4*  75.2  8.5  242.4  59.5  72.8 

Units:  10 (  9>  N/m2 . 

*  The  value  of  57.3  appearing  in  Ref.  (1)  has  been  changed  so 
that  the  relation  c66  =  (ell  -  cl2)/2  holds?  ell  and  c66  are 
directly  measured  and  hence  are  more  accurately  known  than 
cl2. 

TABLE  11.  PIEZOELECTRIC  STRESS  CONSTANTS. 

el5  e22  e31  e33 

3.76  2.43  0.23  1.33 

Units:  C/m2. 
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TABLE  12.  DIELECTRIC  PERMITTIVITIES  AT  CONSTANT  STRAIN. 

(eps)Sll  (eps)S33 

392.  247. 

Units:  10<-12)  F/m. 

OUTPUT  VALUES  FOR  LI  NB  03 

The  input  values  from  Tables  10,  11,  and  12  were  used  to  compute 
the  remaining  elastic,  piezoelectric,  and  dielectric  quantities  for 
lithium  niobate  in  the  manner  discussed  in  prior  sections  of  this 
report.  The  results  are  given  in  Tables  13  to  20. 


TABLE  13.  ELASTIC  STIFFNESSES. 


CE 

CD 

CP 

del  cDE 

del  cPE 

del  cPD 

11 

203.0 

218.3 

218.6 

15.3 

15.6 

0.356 

12 

57.4  * 

42.6 

42.2 

-14.8 

-15.2 

-0.340 

13 

75.2 

76.4 

76.5 

1.24 

1.28 

0.0460 

14 

8.5 

-14.8 

-15.3 

-23.3 

-23.8 

-0.539 

33 

242.4 

249.6 

249.8 

7.16 

7.43 

0.266 

44 

59.5 

95.6 

96.4 

36.1 

36.9 

0.833 

66 

72.8 

87.9 

88.2 

15.1 

15.4 

0.348 

Units:  10<  9)  N/m2 . 


*  The  value  of  57.3  appearing  in  Ref.  (1)  has  been  changed  so 
that  the  relation  c66  =  (ell  -  cl2)/2  holds;  ell  and  c66  are 
directly  measured  and  hence  are  more  accurately  known  than 
cl2. 


TABLE  14.  ELASTIC  COMPLIANCES. 


sE 

sD 

SP 

del  sED 

del  SEP 

del  sDP 

11 

5.83 

5.26 

5.25 

0.574 

0.580 

0.00690 

12 

-1.15 

-0.585 

-0.579 

0.568 

-0.575 

-0.00690 

13 

-1.45 

-1.43 

-1.43 

0.0202 

-0.0209 

-0.000727 

14 

-0.998 

0.905 

0.928 

1.90 

-1.93 

-0.0227 

33 

5.03 

4.88 

4.88 

0.142 

0.147 

0.00512 

44 

17.09 

10.74 

10.67 

6.35 

6.42 

0.0756 

66 

13.97 

11.69 

11.66 

2.28 

2.31 

0.0272 

Units:  lot"12)  m2/N. 
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TABLE  15.  PIEZOELECTRIC  [e],  [h] ,  AND  [a]  VALUES. 


e 

h 

a 

15 

3.76 

9.59 

9.81 

22 

2.43 

6.20 

6.34 

31 

0.23 

0.931 

0.966 

33 

1.33 

5.38 

5.58 

Units:  e: 

C/m2;  h  and  a: 

10<  9)  V/m. 

TABLE  16. 

PIEZOELECTRIC  [d] ,  [g],  AND  [b] 

VALUES . 

d 

g 

b 

15 

69.1 

91.8 

92.9 

22 

20.7 

27.5 

27.9 

31 

-0.854 

-3.36 

-3.48 

33 

6.02 

23.6 

24.5 

Units:  d: 

10<-12>  m/V;  g 

and  b:  10 3) 

m2/C. 

TABLE  17. 

DIELECTRIC  (eps) 

VALUES . 

del  (eps)TS 

11 

392.0 

752.6 

360.6 

33 

247.0 

254.6 

7.61 

Units:  10 

<'12>  F/m. 

del  (eps)TS  =  del  (chi)TS 

TABLE  18. 

DIELECTRIC  (chi) 

VALUES . 

(chi)T 

del  (chi)TS 

11 

383.1 

743.7 

360.6 

33 

238.1 

245.8 

7.61 

Units:  10 

(-!2)  p/in> 

del  (chi)TS  =  del  (eps)TS 
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TABLE  19.  DIELECTRIC  (bet)  VALUES. 


(bet) S  (bet)T  del  (bet)TS 

11  2.551  1.329  -1.222 

33  4.049  3.928  -0.121 

Units:  10 (  m/F. 

TABLE  20.  DIELECTRIC  (zet)  VALUES. 

(zet)S  (zet)T 

11  2.610  1.345 

33  4.199  4.069 

Units:  10 ( * 1 2 3 * * * * *  9)  m/F. 


del  (zet)TS 

-1.265 

-0.130 


CONCLUSIONS 

This  report  provides  formulas  interrelating  the  coefficients  that 
appear  in  the  several  alternative  sets  of  constitutive  equations 
involving  the  elastic,  piezoelectric,  and  dielectric  properties  of 
crystals.  These  are  then  specialized  for  crystals  of  class  3m; 
using  measured  values  reported  for  lithium  niobate,  numerical 
values  of  the  elements  of  the  polarization  matrices  are  calculated. 
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